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Abstract 

We consider the problem of how to construct the curvature perturbation £ to non- 
linear levels, which is expected to evolve time independently on super-horizon scales; 
in particular we concentrate on the situation where the universe is dominated by a 
perfect fluid. We have used a low energy/long wavelength effective action to model the 
fluid sector. Different from previous work, our approach assumes neither the absence of 
vector and tensor perturbations nor "local homogeneity and isotropy" . As a corollary, 
we also show that the nonlinearly defined graviton field jij is conserved outside the 
horizon in the same manner as £ is. 

1 Introduction 

Despite the significant success of the linear perturbation theory in cosmology over the last 
few decades, more and more efforts have been devoted in recent years to understanding the 
cosmological perturbations beyond the linear level — both their origins and evolutions. This 
is partially motivated by the extraordinary improvement in the accuracy of experiments; 
for instance, both the CMB data and the large scale structure data have indicated that 
there might be a detectable non-Gaussian feature in the primordial curvature perturbation. 
Therefore to fully exploit the data, one has to study various cosmological perturbations at 
non-linear level. 

Among these quantities, enormous attention has been paid to the curvature perturbation 
(. Originally it was defined at linear order as the scalar metric perturbation on the uniform- 
density hyper-surfaces by Bardeen [1] and was shown to evolve time-independently outside 
the horizon (cf. for instance |2J). Many studies have been done recently in seek of extending 
the construction of a conserved £ to non-linear level in different contexts. Roughly speaking, 
there are three classes of methods employed to achieve that, which are summarized as follows: 

• The first one is the standard perturbative approach — a perturbative expansion in fields 
|12j . It is straightforward and the equations governing the evolution of ( are valid on 
all scales. However the price to pay is that the equations usually become cumbersome 
once one goes beyond the linear order in fields and there is a lack of systematic ways 
to demonstrate the conservation of ( on super-horizon scales to arbitrary orders in the 
field expansion. 
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• The second class of approaches instead focuses on the spatial gradient expansion - 
i.e. a expansion in powers of a = k/aH, where k denotes the wavenumber, — since 
the super-horizon evolution of ( is usually of most interest. At the leading order, it 
corresponds to the separate universe approach, HD] investigated the cases with an 
inflating background. There the authors showed that in single field inflation models, 
( can be defined as the (only) dynamical scalar metric perturbation in the unitary 
gauge (in which the matter field is demanded unperturbed) and that such a ( was 
conserved to all orders in fields on super-horizon scales. The proofs were based on a 
Lagrangian formalism, for the matter contents in these cases were specified by explicit 
matter actions. On the other hand, the cases with a general FRW background were 
studied in [HHHHS] (known as the SN formalism). It was shown that the conservation 
of the nonlinear £ on large scales was assured if the pressure of the matter content was 
only a function of the energy density or if the matter content could be modeled by a 
single scalar field. As opposed to the aforementioned field theoretic formalism used in 
inflation scenarios, these analyses relied on the field equations of motion; in the context 
of the Einstein gravity the energy conservation alone would be sufficient. 

• A third approach was proposed in [5], E] by invoking the purely geometrical description 
of the curvature perturbation. The covariant quantity they constructed there would 
satisfy an exact, non-perturbative, valid-to-all-scale conservation equation and would 
eventually reduce to usual ( on large scales. The equivalence between this covariant 
formalism and the SN formalism was established in [HI [T3] . 

In this paper, we will present a systematic way of constructing the nonlinear conserved 
( in a universe dominated by a perfect fluid. A low energy effective action for an ordinary 
perfect fluid [U [9] will be employed to model the matter content. So our method here is 
parallel to that used in the single field inflation cases [7].[TU]. On the other hand, although this 
scenario was investigated to some extent in previous literature, we believe that there is still 
some novelties in our approach: i) It does not assume that on a sufficiently large smoothing 
scale, the universe looks locally homogeneous and isotropic (like a FRW patch), as was a 
crucial assumption in (TTJ 0]; in fact, such a feature, though can be derived, does not play 
a vital role in our construction. And ii) our approach works in the presence of vector and 
tensor perturbations. Moreover it enables us to construct the vector and tensor counterparts 
of (, — i.e. the nonlinear vector and tensor perturbation that evolve time-independently 
outside the horizon. Q 

2 Setup 

We begin with a brief review of the low energy effective description of an ordinary perfect 
fluid system in Minkowskian spacetime, which was proposed and developed in [HI El [161 E]- 
For the purpose of this paper, we focus on a fluid with no conserved charges and neglect 

lr The conservation of the vector field does not contradict the usual intuition that vector modes decay in 
the absence of anisotropic stress tensor. As we will see later, due to the symmetry requirements of a fluid, 
any constant (in time) transverse vector field configuration is not physical and hence can be set to vanish, 
leaving only the decaying configuration. 



2 



all the dissipative effects. The readers already familiar with our notations are welcome 
to skip directly to Section 2.2 in which the effective fluid action is used as the matter 
action in cosmological context and in which some formulas are derived in preparation for the 
construction of the conserved (nonlinear) quantities on large scales. 



2.1 Effective Theory for Perfect Fluids in Flat Spacetime 

Our goal is to construct a Lorentz invariant low energy effective theory for an ordinary 
perfect fluid system. For an ordinary fluid system with no conserved charges, we can specify 
as the long wavelength degrees of freedom the comoving coordinates of fluid volume elements, 
parametrized by 7 , with 7 = 1,2, 3. So at any fixed time t, the physical position occupied 
by each volume element is given by x(0 7 ,t). In this description, known as the Euclidean 
description, the physical spatial coordinates x % serve as dynamical fields while t and (f) 1 are 
analogous to world sheet coordinates. 

However, it is often more convenient to use the inverse functions <f) ! (x, t) as dynamical 
degrees of freedom (known as the Lagrangian description), since the spacetime symmetry 
can be straightforwardly implemented — we simply demand (p 1 to transform as scalars under 
Poincare transformations. Also we are allowed to choose the comoving coordinates in such a 
manner that when the fluid system is at rest, in equilibrium and in a homogeneous state at 
some given external pressure, cj) 1 = x 1 — in the field theoretical language, this is equivalent 
to specifying the ground state of our theory to be 

(0 7 > = (1) 

What are the other symmetries, in addition to the Poincare invariance, required to make 
the system behave like an ordinary fluid? Notice that our choice of the ground state (jTl) breaks 
both spatial translational and rotational invariance. In order that the energy momentum 
tensor of the system in equilibrium remains homogeneous and isotropic, as is indeed the case 
for a fluid, we shall impose internal symmetries to compensate the spontaneously broken 
spacetime symmetries. More precisely, we demand that the theory be invariant under such 
internal transformations (i.e. all fields evaluated at the same spacetime point) as: 

% : (j) 1 -»> (j)' 1 = (j) 1 + a 1 , a 7 constant (2) 
TZi : (j) 1 (f)' 1 = 0[j<p J , e 50(3) (3) 

It is easy to show that the background configuration (jl]) is invariant under the diagonal 
translation and rotation, respectively, which are defined as a linear combination of the spatial 
and internal transformations : Td = T s + Ti and IZj, = TZ s + TZi] it is these residual symmetries 
that correspond to the homogeneity and isotropy of our background configuration. 

Moreover an ordinary fluid is insensitive to incompressional deformations — it costs no 
energy to displace fluid volume elements if they are not compressed or dilated. Expressed 
in terms of a symmetry requirement, we demand that the theory be invariant under volume 

2 For space filling fluid, there is a differeomorphism between the comoving coordinates cj) 1 and physical 
spatial coordinates x 1 so that for any fixed time t, the inverse functions 4> 1 {x 1 1) exist. 
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preserving differomorphisms of the comoving coordinates, defined as 

V : 7 -> £ 7 (0), with det ^7 = 1- (4) 

So now we are ready to construct the effective action for an ordinary fluid from 7 's, 
compatible with the symmetry properties mentioned above. It is organized as a derivative 
expansion. Since the internal translation symmetry ^ mandates each field to be accom- 
panied by at least one derivative, at the leading order in the derivative expansion, the only 
invariant is 

B = det(B IJ ), with B IJ = d^d^tf (5) 
and hence the effective action takes the form of 

SVi u i d = Jd'xF(B) (6) 

where F is a generic function, which, as we will see, characterizes the equation of state of 
the fluid in question. 

To illustrate that this effective action indeed describes a perfect fluid, we need to check 
that the (conserved) energy moment tensor takes the famous form: 

Tpv = (p + p)u^u v + p rj^ v (7) 

Before doing that let's work out the four-velocity field u^(x) for the fluid system in our field 
theoretical language, which by definition is given by 

= -^V(x) ee u»{x)d^{x) (8) 

ar 

where r parametrizes the streamline and the derivative with respect to r vanishes because 
the comoving coordinate (or the label) of each fluid volume element is fixed. By solving the 
above equation for w M , we obtain 

1 ^ d^dptfd^* (9) 



where e is the 4d Levi-Civita symbol, with the convention 60123 = — e 0123 = 1. The normal- 
ization and overall sign of are chosen such that u^u^ = — 1 and u° > 0. 
The energy momentum tensor following from the effective action (J6| reads 

= -2F'{B)B{B- 1 ) U drfdvtf + V^F(B) . (10) 

With the aid of Q, the energy momentum tensor above indeed can be recast into the perfect 
fluid form ([7]), if we identify p = —F(B) and p = F(B) — 2F'(B)B. This also justifies our 
preceding claim that the generic function F determines the equation of state for fluids. For 
instance, for an ultra-relativistic fluid with p = p/3, one has F(B) oc B 2 ^ 3 . 

In the rest of this subsection, we will consider small fluctuations about the homogeneous 
equilibrium background configuration (jl|); they are associated with Goldstone excitations. 

(f) 1 = X I + 7C I {x) . (11) 
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Not all the it 1 fields feature propagating wave solutions. This can be seen by expanding the 
effective Lagrangian ^ to the quadratic order: 



£ (2) = \v>o (vri - cl0v L y) + Iwoft (12) 



where til and ttt are the longitudinal (curl-free) and transverse (divergence-free) components 
of 7T 7 : 7 

^ = ^pvr L + 4 (13) 

and where w is the equilibrium enthalpy density p + p. The (squared) speed of sound c 2 is 
given by 

2F"(B)B + F'(B) 



, dp 

c 2 - 



dp 



(14) 

B=B 



B=B F'(B) 

Indeed, we see that only the longitudinal Goldstone field ttl admits the standard propa- 
gating mode at a finite speed c s , with a dispersion relation u = c s k, while the other Goldstone 
field n T has a degenerate dispersion relation uj = and thus does not propagate. For this 
reason, we usually interpret tt l as the sound wave d.o.f. and n T the vortex d.o.f. . 

It is worth pointing out that after spontaneous breaking, the spacetime symmetries get 
mixed with the internal ones, so the three Goldstone fields n 1 transform as a vector field 
under the diagonal 5*0(3). n And after the decomposition (13), n L can be regarded as a 
scalar field and ttt as a transverse vector field. 



2.2 Cosmological Models with Perfect Fluids 

In cosmology, the matter contents are often modeled as perfect fluids. Fortunately, in our 
field theoretic language, we can straightforwardly write the action of the cosmological model 

as 

S = SeH + SVluid (15) 

where the first term on the r.h.s. is just the usual Einstein - Hilbert action for gravity and the 
matter action Spiuid is given by the fluid effective action (|6]), with the flat metric 77 replaced 
by a cosmological spacetime one g and the measure by \f—g d A x . We also assume that 
the metric fluctuates around a flat FRW background g^ v = diag{— 1, a(t) 2 , a(t) 2 , a{t) 2 }. The 
number of dynamical degrees of freedom in question is counted as follows. In addition to the 
three Goldstone fields 7r J — which, as argued in last subsection, are eventually regrouped 
into one longitudinal scalar and one transverse vector under the residual £0(3) group — of 
the fluid sector, the gravity sector introduces one extra 2— polarized d.o.f. - the spin— 2 
graviton, a traceless transverse tensor field. 

Perhaps a more transparent way to see this is by removing the gauge redundancy of the 
gravity sector. Using the ADM variables, one knows that only the spatial part of the metric 
is dynamical, which can be cast into the most general form of 

9ij = = aitf exp (ASy + didjX + diCj + djQ + Aj) (16) 

3 For this reason in the rest of the paper, we will not distinguish the spatial label "i, j, . . . " from the 
internal ones "I, J, ..." . 
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with C i} Dij satisfying diCi = and diD^ = Da = 0. If the four gauge conditions of 
coordinate transformations are chosen to set A = 0, \ = 0, C, = Q — which is known 
as the spatially flat slicing gauge (SFSG), — we are left with til, ttt, and Dij as the 
dynamical degrees of freedom, which characterizes, respectively, the scalar, vector and tensor 
cosmological perturbation. 

For the interest of this paper, it turns out to be most convenient to work in another 
gauge — called the unitary gauge (UG), — in which all the perturbations are absorbed into 
the metric, leaving the matter fields unperturbed: cj) 1 = x 1 . Said differently, the spatial 
coordinates are chosen to coincide with the comoving ones. Meanwhile the temporal gauge 
freedom is used to determine the time slices such that the scalar quantity B (defined in eqn. 
(J5|) remains unperturbed on each time slice: B(t) = a(t)~ 6 . 

The unitary gauge leads to many conceptual and computational simplifications. First 
of all, with our choice of the spatial coordinates, the worldlines of fluid volume elements 
coincide with the threads x % = constant. Indeed, this can be seen by considering the spatial 
components of the velocity field, which vanish since u 1 oc e Iabc d a (f) 1 db<j) 2 d c <p 3 = 0. Second, 
our choice of time slices coincides with the uniform density slices, for the energy density of 
the fluid, given by p = —F(B), is only a function of B and it is a constant on each time 
slice. Using the ADM variables, we can parametrize the metric as 

ds 2 = -Af 2 dt 2 + h i:j (dx l + N l dt) {dx j + N j dt) (17) 

and the inverse metric g^ u as 

^oo 1 „oi io jj j- N' l N J 



where the spatial metric hij takes the form of (16) and K 11 is the inverse spatial metric: 



h hkj = 6j. One can show easily that in the UG B — g and that 

where the index of N 1 (Nj) are lowered (raised) by hjj (h IJ ). Therefore the time slicing 
condition can be expressed as 

B(t) = a(t)- & 3A + V 2 X = log(l-^). (20) 

Now the number of dynamical d.o.f. in the UG can be counted as follows: since M and 
N 1 are just auxiliary fields and can be determined in terms of hy after algebraically solving 



the imposed the constraints, (20) implies that the two scalar functions A and x are n °t 
independent. So henceforth we shall always regard \ as being expressed (perturbatively) in 
terms of other perturbations in the metric and select A, C< and as the dynamical fields, 
the number of which is in agreement with that in the SFSG. 



4 We show in Appendix |a] that this can always be done via a perturbative construction. 
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We are ready to solve in the UG the constraint equations 5S/5N = 0,5S/5N l = and 



the time slicing condition (20) up to the linear order in fields and to obtain the UG quadratic 
action: 



c(2) 

kJrp 



dt 



a{tf 



k 2 

D 2 -—D 2 - 

V a 2 V 



S 



(2) 



-M 2 x Ha(tf 



dt I a(t) 3 : 

k 1 - AHa 2 /k 2 



C 2 



s 



(2) 



dt I a(tf 

k 



VM^Hait) 2 
A{3Ha 2 - k 2 ) 



A 



H 
H 



A 



9M&H 



1 + 



3HH 



A 



(21) 
(22) 

(23) 



d 3 k 



where the symbol is abbreviation for J 
equation 



and where we have used the Friedmann 



-F(B) = 3M^H 2 , p + p=-2F'(B)B 



-2M^H. 



In the long wavelength limit a = k/aH — > 0, the actions 23 (reduce to 



s< 2 > = + s, 



A 2 ) 
v 



dt I M^aity 

k 



4 



A 



(24) 



(25) 



A few comments are in order: 

^2 



i) The (squared) speeds of sound for fluids in curved 
spacetime are & s = dp/dp\§ = — 1 — H/3HH and c\ = 0, respectively, for the longitudinal 
and transverse Goldstone excitations. For c s to be sub-luminal, we need 1 < —H/3HH < 2; 
for fluids with constant w = p/p, this condition corresponds to the scale factor evolving as 



a(t) oc t n } with | < n < |. In particular, this implies that our cosmological model (15) 



with a perfect fluid can not be a consistent inflationary model. Indeed to have inflation, 
we need to relax some symmetry requirements for our system, using less symmetric object 
(e.g. a solid [18]) as the matter content, it) The IR quadratic action (25) implies that we 
should really treat A ~ kCi ~ as of the same order in the spatial gradient expansion. 
Also despite the deceptive appearance, V 2 x is not necessarily of higher order than A in the 
spatial gradient expansion; in fact as we will show soon, they are both of the leading order. 



Thus the scalar, vector and tensor metric perturbation in (16) are all of the same order in 
the spatial gradient expansion, i.e. 



A ~ k 2 x ~ kd 



Di 



(26) 



3 Construct the Conserved Curvature Perturbation to 
Non-linear Order 

As mentioned in the previous section, if the matter content of the universe can be mod- 
eled as an ordinary perfect fluid, there is only one dynamical scalar field. In the UG it 
is parametrized by A(x), the (dynamical) scalar perturbation in the spatial metric (eqn. 



(16)). In this section, we will show that the evolution of A will remain time-independent 



as long as the mode is outside the horizon. The conservation is preserved up to all orders 
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in the field expansion; in particular, at the linear order A coincides with (twice) the usual 
(linear) curvature perturbation. Thus we shall define our linear curvature perturbation ( by 
C = A/2. 

The proof of the conservation of £ on large scales will proceed in several steps as follows: 
Firstly we expand the action (15) to the leading order in the spatial gradient, i.e. to the 



order O(a ), while keeping all orders in fields. We show that the IR action starts with terms 
involving at least two time derivatives, each of which acts on a different field, and hence the 
constant (time-independent) configurations of A,Ci,Dy are permitted as solutions to the 
classical equation of motion. And secondly we show that these constant configurations are 
stable under small fluctuations — i.e. they are indeed attractors. As we will see the proof 
for the scalar, vector and tensor perturbations are essentially identical. 
Step 1: To show the conservation of A,Ci,Dij, we expand the Lagrangian (15) up to the 



first order in the temporal derivative and to the zeroth order in the spatial gradients, while 
keeping all orders in fields. The main challenge of doing this is to express the non-dynamical 
quantities, such as Af, N l and \i in terms of the dynamical ones, by using the time slicing 



condition ( 20 ) and the constraint equations which are given by 



5S _M 2 X 



R (3) + 6H 2 + 2H Tr M + -Tr 



+ ViN j ( 



SS 



d 



dt dt 



2H Tr M 

" a7 



F(B) + 2F'(B)B 
, V k N k 



N*Ni 



- 2F'(B)B 



2M M 

UN, 
N 2 - N k N k 



+ V, 



M 2 
1 



N i Ni 



2M V dt 



+ 



- ( TtmV 



(27) 



2A7 



(2f 



where is the spatial Ricci scalar constructed from hij, V the covariant derivative com- 
patible with h^, and the matrix Mjj = \og(hija~ 2 ) . The ". . . " in the first equation stands 
for terms that are regular in the limit a — > — the form of which is irrelevant in the analysis. 



Inspecting the constraint equation (28), one finds that N l starts at least at the order 0(a). 



Thus, as long as the super-horizon modes are concerned, all the terms in (27) with N l or 



V as well as RW) can be neglected, which leads to a crucial fact that SAf = M — 1 in the 
super-horizon regime starts with terms involving two time derivatives: 



8N 



1 Tr ( — P - M —p M 
2AH 2 [dt dt 



-,1/2 



a 



(29) 



Moreover, we see that for each field, there is at most one time derivative acting on it. 



On the other hand, the time slicing condition (20), in this long wavelength limit, reduces 



to 



Tr M = 3A + V 2 x = O (a 2 ) 



(30) 



which verifies our claim that V 2 x is of the same order as A in the spatial gradient expansion. 
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It then follows immediately that in the super-horizon regime, the Lagrangian (15) be- 
comes 

lim C ~ 2a{tfNF{B) ~ -Qa{tf M^H 2 + V C n {<p a ) (31) 

k/aH^O J 

n>2 

In the last step we have denoted the dynamical fields — A,Ci, Dij — collectively by <ft a and 
the number of fields contained in C n by the subscript "n". Notice that schematically £ n 's 
take the form of 

C 2 ~ Q2{t)Gf(j) a fa, C n>2 ~ Q n {t)Gt M {<i>)<i>a fa---Mi (32) 

with Q n (t) being a function of time consisting of a(t),H etc.. Besides the irrelevant field- 
independent term, this long-wavelength Lagrangian starts at two time derivative level and 
hence the e.o.m. following from it reads 

f x {<j>, 4>)fa + f 2 (<f>,4>) b J b = (33) 

Therefore it indeed admits A, C{, = constant as solutions to the classical equation 
of motion. Since this IR Lagrangian contains all orders in fields, the conservation of these 
fields on large scales must be preserved nonlinearly. 

Step 2: Now we show the solutions A, C;, = constant are actually attractors. We just 
work on the scalar case, since the analysis for the other two is identical. Writing A as 
A = Aq + SA, owing to the constancy of Aq, the quadratic action for the fluctuation SA in 
the long-wavelength limit takes the same form as that for A, which is given by 

6S<?> = j 'dtJ^a(t) 3 6A 2 

from which follows the linearized equation of motion for SA 

5A + 3H5A = Q (34) 

It admits two general solutions — one decaying mode and one constant mode: 

f dt 

b~A\ = / , 5A 2 = const. (35) 

J a(ty 

This thus confirms that the solution A = constant is an attractor. 

As we saw, our approach is more powerful in some aspect than those in previous literature, 
for it enables us to show that the vector perturbation Cj and the tensor perturbation 
are also conserved on super-horizon scales in the same manner as their scalar counterpart 
A. And in next section, we will show that the nonlinear, conserved ( = A/2 and we 
constructed here agrees with that in jH [12] , up to a global spatial coordinates redefinition. 

4 Discussion 

We have constructed the curvature perturbation ( (to nonlinear level) via the following steps: 
i) choose a coordinate system such that the spatial coordinates comove with the fluid and 



9 



that equal time slice coincides with the uniform density slice. And ii) define ( to be (half) 
the coefficient of the term in log(a~ 2 gij) that is proportional to 8{j. That is, 

C = ^V" 2 (c% V 2 - dtdj) log [a- 2 9ij ] (36) 

We showed that the ( such defined evolves time-independently in the super-horizon regime 
to all orders in fields. 

Our result of the conservation of the non-linear curvature perturbation ( on large scales 
agrees with previous literatures jH [T2] . But in our analysis we did not neglect the presence 
of the vector and tensor perturbations, nor did we use the property of local homogeneity 
and isotropyj^] As we showed explicitly, the spatial metric in UG on large scales (er — y 0) is 
given by 

1 



gij = a(t) exp 



6( ( — hkj ) + vector + tensor 



(37) 



At a glimpse, one may think that even though the vector and tensor perturbations are ne- 
glected, this metric is anisotropic, which is inconsistent with our intuition about cosmological 
fluids that no anisotropy in the metric will be generated if the fluid is free of anisotropic 
stress. However, this puzzle results from our non-conventional choice of coordinates: we can 
perform a further coordinate transformation (from the UG) 

t -> t, x 1 y 1 = x 1 + £*(x) (38) 

such that 

<p rI (t, y) = y 1 + ^(y), and 9ij (t, y) = a(i) 2 e 2 ^ exp ^(y), with d a ' tJ = Ju = (39) 

where the time-independence of £ 7 (af) follows directly from the time-independence of A,Ci, 
and Dij and where the function £ /7 (y) is obtained by inverting the 3-diffeomorphism y^x) 



x 1 + £ I (x). Notice that the coordinate transformation like (38) alters neither the time slices 



nor the threading of the spatial coordinate, the latter of which is because 

U ^ JaB-ya i/la i/2a j/3 _ n „,i _ „ „,0 



u oc e 



^d^d^d^' 3 = 0, u 1 ! = g 0I u° = h^Wu ~ 0{a). (40) 



Therefore this implies that in the new coordinate system and in the super-horizon regime, 
both the spatial metric and (c.f. eqn. (|7|) — which is reduced to = pg^ on large 
scales — are isotropic, were the tensor perturbation ignored, and homogeneous since in the 
limit cr — y the spatial dependent of £' and 7' is negligible. 

That is, although convenient for computational purposes, the UG coordinate is not capa- 
ble of exhibiting the local homogeneity and isotropy in the long wavelength limit. Fortunately 
there exists a new coordinate system, which is indistinguishable from the UG coordinate by 
only inspecting the physical quantities such as p,p,u l etc. and in which the existence of 



5 In fact having fixed the UG by (f) 1 = x 1 and B = a(t)~ 6 , we have depleted all guage freedoms for 
choosing a coordinate system. Furthermore, as we argued earlier, we should not think the terms in with 
"superficially" more derivatives such as didjX, diCj are of higher order in the spatial gradient expansion 
than ASij. As a consequence of which, as long as the UG is chosen, we are not entitled to claim any property 
about on large scales. 
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local homogeneous and isotropic patches on large scales can be verified (rather than being 



assumed, as were in previous literatures). The (' in (39) is also time-independent on large 
scales and it agrees perfectly with the definition in [U IT2] . 

As we showed in last section, on super-horizon scales, the vector field Ci evolves time- 
independently. However, this does not correspond to a physical vector mode, as can be 



understood as follows. Transform to the coordinate system specified by (39). Owing to 
the internal volume-preserving-diffeomorphism invariance Q of the fluid system, the con- 
figuration in eqn (39) is in fact symmetric to the one free of vortex degrees of freedom 
(divergence-free vector modes), either in the matter fields or in the dynamical metric com- 
ponents. That is, as the leading contribution, the constant (in time) Ci configuration has no 
physical significance; the well known decaying vector modes in the universe dominated by 
perfect fluids come from subdominant, time-dependent terms in in the UG. 

It is known that a perfect fluid free of vortex degrees of freedom (divergence-free vector 
modes) has an equivalent effective description involving only one scalar field - - i.e. the 
P(X) = P((dijj) 2 ) description [3 [191 120] ■ Then we can apply the method in Ref. [3 [TO] 
directly to a non-inflating background and define the curvature perturbation ( as the (only) 
dynamical scalar metric perturbation in unitary gauge: 

ip(t,x) = ip (t), g i:j = a 2 e 2C exp{7 ii } (41) 

Details of the proof of the conservation of ( on super-horizon scales in the P((dip) 2 ) context 
are collected in Appendix IB] The nonlinear ( constructed via this logic is identical to the £' 



field in (39 ), hence agrees with our definition up to a global reshuffle of the spatial coordinates. 



5 Conclusions 

In this paper we consider the problem of how to extend the definition of curvature perturba- 
tion ( to nonlinear levels; in particular we concentrate on the situation where the universe 
is dominated by a perfect fluid. We have used a low energy/long wavelength effective ac- 
tion to model the fluid sector and constructed the nonlinear ( as follows: i) Fix the gauge 
such that the comoving coordinates of the fluid coincide with the spatial ones and that the 
constant time slices are the uniform density slices; and ii) define ( = A/2 where A is the 



isotropic scalar part in the spatial metric (16). We have shown this £ (as well as its tensor 



counterpart) is conserved outside the horizon to all orders in fields. 

Although this topic has been investigated to some extent in previous literature, we believe 
our approach is novel and has its own merits. For instance, in discussing the conservation 
on super-horizon scales of ( (the scalar mode), we do not need to neglect the vector or 
tensor modes; on the contrary, we illustrate that there are vector (non-physical) and tensor 
counterparts of ( which are conserved in the same manner as ( is. Moreover, our proof does 
not assume/rely on that the universe looks locally like FRW patches (local homogeneity and 
isotropy) on sufficiently large scales, which was the key input in previous work. Nevertheless, 
we have shown our definition of ( agrees with that in Lyth et al. (2005) and Malik et al. (2004) 
up to a global reshuffle of the spatial coordinates . 
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Appendix 

A the Spatially Flat Slicing Gauge (SFSG) to All Or- 
ders in Fields 

In this section, we show that starting from an arbitrary gij(x) = + 8g^, we can always 
achieve the non-linear SFSG, specified by 

gij(x) = a(t) 2 exp7 ii (x), with 7^ = <%7y = (42) 

via some appropriate gauge transformation x^ — > y^ = x^ Here we only concentrate 

on the most relevant case to our analysis, in which g^ u = diag{ — 1, a(t) 2 , a(t) 2 , a(t) 2 } and 
5g is treated as small fluctuations around the FRW background. We prove this by using a 
perturbative construction. 

First let us invert the function y^(x) and denote it as 

00 

x»{y) =y» + «"(y), with v»(y) = £ vf n) (y) (43) 

n=l 

where ui\ is assumed to be of order (5g) n . Under this gauge transformation, the spatial 
metric transforms as 

9a{x) -> gij (y) = yy dyi J yy dyj J g afl (y + v) (44) 



M« = a - 2 S 9ij + divL + djvfa + 2Hv° (l) 5 l3 (45) 



Consider My = log (a 2 gij). At the first order in Sg, it is straightforward to work out M\- (y), 
which is given by 

\r< u ,,- 2 A,, ;)., i 

Notice that any tensor function 5gij(x) = 1, 2, 3) can always be put in the form of 

5 gij = a{tf (ASij + didjB + diCj + dfi t + D y ) (46) 
with d and T>ij satisfying dfii = T>a = diVij — 0. [j Thus by choosing 

v° {1) = -A/2H, tf 1) = -B/2, 4 (1) = -Cj (49) 



6 To see this, we define A and B to be the solutions of 



3A + V 2 B = a~ 2 Sg u and V 2 A + V 2 V 2 B = a- 2 d i d j 5g ij , (47) 



then define Ci as the solution of 
and then use eqn. ( 46 ) to determine 2?,j 



9j (.4 + V 2 B) + V 2 Cj = a- 2 d l 5g ij (48) 
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where v s and v l T are the longitudinal (curl-free) and transverse (divergence-free) components 
of Vf vl^ = diV?s + fy( n ), we find that M-p = T>^ — i.e. only the transverse traceless part 
survives. 

~ (2) 

At the second order the expression for Mh is quite cumbersome, but that doesn't bother 
us much. Since it takes the form of 

M^f = Sij + div\ 2) + d 3 v\ 2) + 2Hvl 2) 5ij, (50) 

by the same argument as above, we show that all the fields in can be set to zero except 
the TT one. And this process can be repeated to all orders in 5g. 

Now let jij = Yl^Li \ which apparently is also transverse and traceless. The spatial 



metric in new coordinate system indeed takes the form of eqn. (42). 



B Nonlinear Conserved ( in P(X) Fluids 

A perfect fluid free of vortex degrees of freedom has a dual low energy effective description 
involving only one scalar field if). In the cosmo logical context, the action is given by 



S = S E h + J v^gP(X), X = -djl>d v il>gr ■ (51) 

The energy momentum tensor can be obtained: 

Tfa, = 2P'{X)d^d„if) + P(X)g^ (52) 



Noticing that = d^ipX 2, can be recast into the standard perfect fluid form (|7) if we 
identify 

p = 2XP' - P , p = P . (53) 

Assuming that the background configuration of the scalar field is only time dependent 
(if)(x)) = ipo(t) and that of the metric takes the usual FRW form, the Einstein equations and 
the scalar field equation of motion for unperturbed configurations read 

3M^H 2 = -P + 2P'ipl (54) 
M^H = -P'ipl (55) 

(> + 2P'Vo) % + SHP'i) = (56) 

Now let us fix the gauge. For the sake of defining a nonlinear (, the unitary gauge is of most 
convenience: 

if)(x) = if) (t) , gij = a(t) 2 e 2C exp7 ij , (57) 

We shall show the £ such defined is conserved outside the horizon. We adopt similar logic 
as before: first show that £ =constant is a solution in the long wavelength limit, and then 
show that this solution is actually an attractor. 



Using the ADM variables, the effective action (51) becomes 



S = jNVh (MflOO + *& {E j El - E 2 ) + P(X)) (58) 
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We can then obtain the constraint equations by varying this action with respect to Af and 



Af 2 
1. 



6(H + C) 2 + 4(H + C)V 4 iV* - -Tr 



d 

dt' 



e 7 — e 7 ] + -V*iV,- V 3 Nj 



(59) 



+ -VjN 1 V J Ni - ViN 1 VjN 1 - 



+ 




4 4 pl 



Mpj Af 2 



+ V 




-^r- r ViN j - -^r- f V j Ni 

2Af 2Af 



(60) 



From equation (60), we can see that ViN 3 ~ C (or 7), i.e. ViN J ' is of the zeroth order in the 
spatial gradient expansion (a = k/aH). Thus solving SAf = Af — 1 via equation (59) up to 
order O(a ), we have 

-1 



SAf 



Mp Y ( P — P'ffl + 2P> 4 ) ( ZH( - HViN* ) + 0(C, Ci). 



(61) 



Plugging this into the action (58) and expanding it to O(a ) while keeping all orders in 
fields, we have 



S = -2A*f j J Vh{3H 2 + H+ (3H( - HV l N i ^j + 0(( 2 , 7 2 , (7) } 



-2M 



pi 



dt 



(aV^)+0(C 2 , 7 2 ,C7) 



(62) 



where in the first equality we have used the background Einstein equations (54) (55) and 
the background scalar field equation (56) and in the second equality we have neglected the 
contribution from ViN 1 since it is a boundary term. Therefore we conclude that the effective 
action in the long wavelength limit starts with terms involving two times derivatives (acting 
on different fields) and, by the same argument as before, that £ = , 7 = are solutions to 
the classical e.o.m.'s. 

Now we show that the solutions ( = constant and 7 =constant are attractors. As before 
we consider the fluctuations around these classical solutions and work out the quadratic 
action for S( and 8j. Owing to the constancy of the classical solutions, the quadratic action 
for fluctuations take the same form as that for ( and 7 themselves, which can be obtained by 
solving the constraint equations (59) (60) up to the linear order in fields, plugging into (58) 
the linearized solutions SAf%, N{ and expanding the effective action to the quadratic order. 
Thus we have 



5S ~ M? 



pi 



a 



£ K 2 - -AdjcY 



+ 



where e 



_c; <r 

-H/H 2 and c s is the speed of sound defined as 
2 _ d P 



cr 



dp 



x 



-41 P'^ 2 ) + 2^ 2 P"^ 2 ) 



(63) 



(64) 
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Thus in the long wavelength limit a — > 0, the linearized equation for 5( (£7) possesses two 
general solutions, one being time independent and the other decaying, which confirms our 
claim that ( = constant and 7 = constant are attractors. 

Before ending this section, let us remark that the nonlinear ( constructed in P(X) context 
agrees with that in F(B) context. Notice that in unitary gauge (57), X = ip (t) 2 (l + 8Af)~ 2 . 
Since on super-horizon scales 8M — > due to the constancy of ( and 7, the constant time 
slices in the unitary gauge coincide with the uniform density slices (for p is a function of 
X only). Similarly, since N l vanishes on large scales for the same reason, we have u l oc 
gOi oc N l — )• 0, i.e. the threading of the spatial coordinates is chosen such that the threads x l 
=constant coincide with the integral curves of the 4- velocity w M (the comoving worldlines). 
Therefore ( defined in P{X) context is identical to (' in (39). 
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